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The aspect of the quasiparticle interaction of a local Fermi liquid, the impurity version of 
f^-based heavy fermions, is studied by the Wilson numerical renormalization group method. 
In particular, the case of the P-singlet crystalline-electric-field ground state is investigated 
assuming the case of UPt3 with the hexagonal symmetry. It is found that the interorbital in- 
teraction becomes larger than the intraorbital one in contrast to the case of the bare Coulomb 
interaction for the parameters relevant to UPt3. This result offers us a basis to construct a 
microscopic theory of the superconductivity of UPt3 where the interorbital interactions are 
expected to play important roles. 
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Heavy fermions containing two f-electrons at each site (f 2 ) have attracted much attention 
for a decade or so. Although uranium compounds such as URu2Si2, UBei3, and UPt3, and 
Pr-based filled skutterudites exhibit rich phenomena such as superconductivity, anomalous 
magnetism, non- Fermi liquid (NFL), metal-insulator transition and so on, a microscopic the- 
ory for these systems does not seem to have been well developed. Coming to think of low 
temperature physics of those f 2 systems, particularly superconductivity, we need to know on 
the basis of a microscopic approach how quasiparticles of f^-based compounds interact with 
each other. In this paper, we employ the numerical renormalization group (NRG) method 1 ' 2 
at the level of impurity problems considering the case of UPt3 3 with the hexagonal symmetry, 
because the local part of quasiparticle interaction is expected to play a crucial role also in the 
lattice problem. Our basic assumption is that the ground state of a crystalline-electric-field 
(CEF) is the f^-singlet. The Pt Knight shift experiment 4, 5 and the magnetic succestibility of 
UPd2Ai3, 6 which has the same crystal structure as UPt3, indicate this assumption is reason- 
able. 

In this paper, we use a two orbital Anderson model as an appropriate one for a system 
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Table I. Irreducible decomposition of direct products under D S h ■ Only the different orbital combina- 
tions arc listed. j z = ±5/2(r 8 ), j z = ±l/2(r 7 ), jz = ±3/2(T 9 ) 

with the f^-singlet ground state under the hexagonal symmetry. The validity of this choice is 
discussed in detail in the following. 

The ^-states split into three doublets in the J =5/2 manifold. When we assume a strong 
spin-orbit interaction, we can neglect the 3=7/2 states and use a j-j coupling scheme to obtain 
the f 2 -states. On the other hand, if the Hund-rule coupling is stronger then spin-orbit coupling, 
the f^-states are constructed in Russell- Saunders coupling scheme. The states with 3=4 split 
into five different symmetry states. 

We assume the T 4 singlet ([|3> - | - 3)]/y/2) is a ground state as mentioned above. 
The problem is how we connect these Hilbert spaces, which are in the opposite limit. When 
we construct the f 2 -states in the j-j coupling scheme using only J=5/2 states, we obtain 
no singlets, because the states in the J =5/2 manifold under the hexagonal symmetry are 
the eigen states of the z-component of the total angular momentum operator J\ ot (matrix 
element of Og is zero for J = 5/2 states 7,8 ). Thus, the direct j-j coupling scheme using only 
the 3=5/2 manifold is not appropriate in this case. In this paper, we use the J=7/2 states 
as intermediate states effectively to obtain f^-singlet states. This stems from the fact that in 
the 3=1/2 manifold, the matrix element of 0% is not zero. This fact generates the interaction 
term, which do not conserve J* ot in the 3=5/2 manifold as we eliminate the 3=7/2 states 
from the consideration of Hilbert space(see Fig. 1). This enables us to obtain a singlet ground 
state and we write it as 

Uf 1= * =1 ftAf = *f=i+li.c. . (1) 

2222 22 2 2 

In eq. (1), we ignore retardation and fj z is the destruction operator of the f-electron in the 

3=5/2 manifold. The magnitude of Uf\ _ 5 _x is expected to be small. A detailed analysis 

22TT 

of this term will be published elsewhere. Other interactions such as the Coulomb interaction 
between the f-electrons and the CEF mixing are renormalized. However, their symmetry could 
not be changed. Note that eq. (1) is the only allowed new interaction under the hexagonal 
symmetry in the effective 3 = 5/2 manifold. 

Next, we restrict ourselves to the J z =±5/2, ±1/2 states of 3=5/2 manifold, because 3 Z = 
±3/2 states have no contributions to the f 2 -singlet states(see Table I). Other states are not 
important in the discussion of the low-energy quasiparticle properties under our assumption; 
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Fig. 1. Simplest diagram that generates the interaction of eq. (1). (J,J Z ) denote the angular mo- 
mentum J and its z-component J z . The straight lines represent the J = 5/2 states and the wavy 
lines correspond to the J = 7/2 states. The dashed line and open dots represent the Coulomb 
interaction and the CEF term, respectively. 
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Fig. 2. f 2 level scheme. Ai and A2 are related to K z and K± in eq. (8) as K z — A2 — 2Ai and 
K± = — A2, respectively. 



f^-singlet ground state. Indeed, many band calculations suggest that the J = 5/2 states are 
dominant around Fermi surfaces for UPt3. 9 ~ n 

Before writing our Hamiltonian, we transform representations from j z to pseudo-spin. This 
is performed using 



/s -> /at, /_5 /A|j (2) 

— ► /fit' A — *■ fsi, (3) 
^ 2 

where fj x and / a „ are f-electron field operators with the z-component of angular momentum 
j z , and a and /i are the orbital and pseudo-spin indices, respectively. In this representation, 
the f 2 -atomic energy spectra are given by an anisotropic antiferromagnetic Hund's-rule cou- 
pling 14 (See Fig. 2). The splitting of A 2 originats from eq. (1) (redefined as Uf[ _ 5 _! = K±/2). 

2 2~~T 

Then our Hamiltonian is 



H 



H c + Hi oc + i^hyb + -^Hund, 
€ kc\. afl Ckafi-, 



(4) 
(5) 



kafj, 



loc 



e fJU<*» + 2 v <* n u n u 

fia a 
+UABj2 n fA» n fB„i 



(6) 
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#hyb = y^^VkgCkanfan + h.C.], (7) 

kfiCt 

-f^Hund = -K Z S Z A -S Z B ^[S^Sg + S^Sg], (8) 



where S^,S^ are the z- and transverse components of the f-electron pseudospin operators 
of the a(=A,B) orbital. The other notations are conventional. In eq. (7) we have assumed 
that each f-electron mixes with the conduction electrons with the same point group symmetry 
and two conduction electron bands (from — D to D) do not mix with each other. Note that 
the exchange interactions among different f-orbitals at the impurity site are not important 
for discussing the f 2 -singlet ground-state model, because f 2 -states include such terms in their 
energy levels. The Hamiltonian (eq. (4)) can be regarded as a "simplified" two-impurity An- 
derson model, although it contains an anisotropic spin-spin interaction and two nonequivalent 
impurities. The word "simplified" means that the conduction electrons belonging to one band 
do not mix to those of the other band after impurity scattering. 

In order to solve the Hamiltonian eq. (4) by the NRG method, we discretize conduction 
electron bands logarithmically and transform eq. (5) into one-dimensional semi-infinite chains, 
following Wilson. 1 Then we write eqs. (5) and (7) as 

oo 

H * = J2I2 A ~ n/2 OfUfn + la»+h.C.}, (9) 
aft n=0 

-^hyb = ^V a [fi la ^f 0af , + foapf-latt], (10) 

otp, 

t a = q + A-^q-A-"- 1 ) 

2y/(l- A- 2 «- 1 )(l-A- 2 ™- 3 )' 

vl = £l*U 2 ' (12) 

k 

where f-i a u = /au ; fna^ ( n > 0) is related to the conduction electron field operator via 
unitary transformations, and we take D as the unit of energy. In this paper, we take the 
logarithmic discretization parameter as A = 2.5. 

The Hamiltonian eq. (4) describes the competition between the Kondo-Yosida singlet and 
the f 2 -CEF singlet. There are two stable Fermi liquid (FL) fixed points: the Kondo-Yosida- 
singlet fixed point and the CEF-singlet one. Between these two FL fixed points, there exists an 
unstable non-Fermi liquid (NFL) fixed point. 14 Around the above two stable FL fixed points, 
we can write the effective Hamiltonian as 

= h*+ [tad la + u a d 2a }+u AB d 3 

a=A,B 

+K z d 4 + K ± 5 , (13) 

where t a ,U a , Uab,K z and K± are parameters that should be determined by fitting the NRG 
spectrum near the FL fixed point whose Hamiltonian is H* (in this case, free Hamiltonian), 
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TV means the number of conduction electron shells in the NRG procedure and O's are defined 

as 

Ola = A^/oa^K/o^ + ^/L^+h-C, (14) 

d 2a = A^a^-l) 2 , (15) 
3 = A^aoVo^-^aoB./oB.-l), (16) 

°4 = — (/oat /OAT - /oAi/oAl) 

x(/ t J BT^osT ~ /obt/obi)> (17) 

A 2 1 . 
°5 = — ^— (/oat/oat/obt/obt 

+/o t Al/oAT/o t BT^' OB i)' ( 18 ) 

where u Q is the potential scattering amplitude of orbital a at the FL fixed point and we use 
the Einstein contraction notation for \i and v. This type of effective Hamiltonian was discussed 
in terms of the two-impurity Kondo problem except for the particle-hole and spin rotation 
symmetry. 12 Using eq. (13) we can extract quasiparticle interactions as 

rpaa _ 2-paaaa (1 ~l~ A ) U a (~\Q\ 

rpAB _ „ „ \ ^ r ABAB _ 1 + A 1 Uab (e) „, 
I< s = z A z B 2^ L WT* ~ 4 f^' ( 20 ) 



tthAB - ^ „ r ABAB 1 + A 1 K z /A 



(21) 



rpAB _ „ „ -pABAB 1 + A 1 if±/2 

f cx = ^r TUT = 4 ^ • (22) 

Here, r"J" 2 2 " 3 3 " 4 is a full vertex part whose energy variables are set to zero (on the Fermi 
surface), z a oc t^ 1 is a renormalization factor of orbital a, and f (J,) corresponds to <r = 
+1(— 1). These results are obtained by comparing the form of pseudospin susceptibility and 
specific heat between NRG and local Fermi liquid theory. Note here that eq. (22) is not an 
exact relation, because the transverse pseudospin components are not conserved (we take the 
z-axis as the quantization axis), so that the contribution from incoherent parts remains finite. 
Although we cannot write F ex only in terms of quasiparticles, in the NRG calculation, we 
think eq. (22) is a reasonable expression. 15 Noted that the interaction F^ B stems from the 
CEF and does not conserve the z-component of the total angular momentum. 

The results of quasiparticle interactions and Wilson ratio are shown in Figs. 3 and 4, 
respectively. Although we have performed calculations for various parameters, we have only 
shown the results for a specific parameter set because other sets also show qualitatively the 
same behavior in the Kondo regime. Note that the level of excited Ts coincides with that of 
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Fig. 3. Quasiparticle interactions vs The parameters used in the calculation are 6a — —0.6, c_b = 
-0.5, U A = U B = 2.0, Uab = 0.4, v A = v B = 0.25 and K± = -0.04. The NFL fixed point is 
K c z ~ —0.02 for these parameters. Note that Ai is zero at K z = —K± = 0.04. 



Fig. 4. Wilson ratio vs K z . R a is the Wilson ratio of orbital a without interorbital contributions. R is 
the total Wilson ratio. The parameters used in calculation are the same as in the Fig. 3. These are 
calculated by peudo-spin representation, so that the results are different from real Wilson ratio. 

T4 at K z = —K± and the level splitting between T§ and T4 increases as K z decreases. The 
interesting point in Fig. 3 is that the interorbital interactions F^ B and F^ B are enhanced 
around the intermediate coupling region (K z ~ 0.01) and F^ B becomes considerably larger 
than the intraorbital interactions. For a large K z , this system approaches the fixed point of 
the anisotropic S = 1 two-channel Kondo model, where the f-electron renormalization factors 
z a (a = A, B) become small. At the NFL fixed point (K z = —0.02), z a vanishes, which means 
the quasiparticle description breaks down. 

In Fig. 4, both Ra and Rb are almost constant and close to 2 (see the caption of Fig. 4 
for the definition of R a ). This means intraorbital quasiparticle interaction is proportional to 
the renormalization factor (F" a ~ z a ~ Tk , where Tx a is the Kondo temperature of orbital 



3.5 - NFL ; 
.o 3.0- i ! 

2 2.5- b 





6/9 



J. Phys. Soc. Jpn. 



Letter 




Fig. 5. f- and conduction electron configuration for the case of (a) T5 ground state and (b) 
ground state. The circles with thick and thin lines represent f-electrons and conduction electrons, 
respectively. The r$ doublet state is constructed by f\^fg^\0) and /Jn/gJO). For the I\ singlet, 

7f [/JiT^si ~ /aj./bt1I ^ 



a). This is a well-known in the two-impurity Kondo problem. 12,13 

There exist four parameters describing low-energy physics, i.e., Tk u (a = A,B), F^ B and 
F^ B . F AB is almost zero in the Kondo regime where the charge susceptibility and orbital 
susceptibility are suppressed. 

A tendency of interorbital couplings, FAB. > reflects the property of the ground-state con- 
figuration of f-electrons (f 2 - F±, T^). If the l^-CEF ground state is the T§ doublet (see Fig. 5(a)) 
and the antiferromagnetic coupling between the f-electrons and conduction electrons (Kondo 
type exchange) is taken into account, we obtain effective ferromagnetic coupling between the 
different orbital quasiparticles, consistent with F^ B < in Fig. 3. The same argument as that 
in the case of the T§ ground state is possible for the ^-singlet ground-state case, see Fig. 5(b). 
However, the region where the Ts-doublet is the ground state, is not valid in the present model 
because of Hund's rule. 16 In the case of f 2 -r 4-singlet ground states, we need not worry about 
Hund's rule, because our f 2 -singlet states have only J = 4 components. In addition, due to the 
reduction of z a toward the NFL fixed point the interorbital couplings F AB ^ are suppressed 
as shown in Fig. 3. Around the intermediate region, the interorbital quasiparticle interactions 
are well enhanced compared with the intraorbital ones (F AA , F BB ), which decrease 
monotonically toward the NFL fixed point. It is interesting that although the overall energy 
scale becomes small as we approach the NFL fixed point, which is a direct consequence of 
z a — > 0, the details of interactions i.e. anisotropies vary. 

The above results imply that the interactions between quasiparticles, in a system with 
multiorbitals (specifically f 2 -singlet ground state), are quite different from a simple model, 
such as a simple-minded multiorbital Hubbard model, in general. If we use the renormalized 
interactions obtained here to tackle the f 2 heavy fermion superconductors, there exists a pos- 
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sibility of obtaining spin triplet superconductivity because of large interorbital interactions. 17 
Several theories of triplet superconductivity induced by Hund's rule have been proposed thus 
far. 18 ' 19 Our results are, however, different from those in the point that the relevant interac- 
tions are those for quasiparticles (in impurity version) and their interactions originate from 
the low-energy CEF level f-electron. Although the contributions of J z = ±3/2 have not been 
included in the present model, we believe they would not change the quantitative features of 
our results. The actual calculation of the superconducting transition temperature of UPt3 is 
very complicated, because of two inequivalent sites in the unit cell and the orbital degrees of 
freedom per site. This will be discussed elsewhere. 

In summary, we have discussed the quasiparticle interactions of a local Fermi liquid with 
the f^-singlet ground state under the hexagonal symmetry. We have found by the NRG method 
that interorbital interactions are much larger than intraorbital ones in the wide range of the 
CEF. This result is in contrast to that in the case of usual Coulomb interactions, in which 
intra-orbital couplings are thought to be much larger than interorbital ones. This may be a 
good starting point to discuss the real lattice systems with two f-electrons per site. 
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